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Abstract. We report on some recent work of the authors showing the rela- 
tions between singular (point) perturbation of the Laplacian and the dynamical 
system describing a charged point particle interacting with the self-generated 
radiation field (the Maxwell— Lorentz system) in the dipole approximation. We 
show that in the limit of a point particle, the dynamics of the system is de- 
scribed by an abstract wave equation containing a selfadjoint operator H m 
of the class of point interactions; the classical Abraham-Lorentz-Dirac third 
order equation, or better its integrated second order version, emerges as the 
evolution equation of the singular part of the field and is related to the bound- 
ary conditions entering in the definition of the operator domain of H m . We 
also give the Hamiltonian structure of the limit model and, in the case of no 
external force, we study the reduced dynamics on the linear stable manifold. 



Dedicated to Sergio Albeverio 



1. Introduction 

Over a century after the discover of the electron by Thompson and the first 
theoretical studies by Lorentz, a satisfying mathematical description of the inter- 
action of point charged particles and electromagnetic field is still lacking. As it is 
well known, the Maxwell-Lorentz system, which correctly describes the bulk mat- 
ter, looses its meaning in the case of a point particle, due to the elementary fact 
that the solutions of the field equations are not regular enough to be evaluated at 
a single point, as required by the particle equation. This classical singularity yields 
directly to the need of infinite mass renormalization and to the difficulties which 
plague the classical theory of the electron. This theory, as emerges in the work of 
Lorentz, Abraham and Dirac among the others, leads to a reduced equation for the 
particle alone, the so called Abraham Lorentz Dirac equation (ALD for short), a 
third order differential equation which embodies the interaction between particle 
and field in an inertial term (mass renormalization) and a third order term, de- 
scribing radiation reaction. The mathematically dubious procedures involved in its 
deduction from Maxwell-Lorentz system, and the presence of unphysical solutions, 
the so called runaway solutions, surround the entire subject of a legitimate suspi- 
cion, and the usual way out is to invoke quantum electrodynamics for the solution of 
the problem. But, in turn, mathematical foundations of quantum electrodynamics 
are far from being clear, and in the opinion of the present authors and many others 
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a solution to the classical problem is a first necessary step toward a clarification 
of the quantum one. In the following (see [2,3] for complete proofs, more details 
and references) we describe the main results obtained in the last years by us in the 
so called dipole approximation of classical electrodynamics, a well known model 
which in essence amounts to linearize the interaction (see the Hamiltonian (3.1)) 
in the Maxwell-Lorentz system. For this model we show that a limit dynamical 
system describing the coupled dynamics of both particle and field indeed exists 
(theorem 3.3), and is an abstract wave equation in which the operator part is a 
point interaction, the elements of the domain of which have a well precise relation 
(see the system (3.5)) with the particle velocity. To simplify the presentation we 
introduce first the special case of a free particle and zero total momentum, where 
the standard delta potential appears, and then we give the general case in which 
an external force is present. To treat this second case a slight generalization of the 
point interaction is needed, where the linear boundary condition of standard delta 
interactions are replaced by an affinc one in which the particle momentum enters 
(see definition 3.2). Such a boundary condition allows to interpret the solutions 
of the limit system as the solutions of the system in which a standard wave equa- 
tion with a delta source is coupled with an ordinary second order equation. This 
equation is a low order version of the ALD equation. 

Our description allows also an immediate Hamiltonian formulation of the limit 
dynamics, a result which solves an old problem in classical electron theory. 

In the last paragraph we study, in the free case, the reduced dynamics on the 
stable linear manifold, also called, in the classical literature, the "non runaway 
dynamics" . In particular we show that the flow of the reduced Hamiltonian system 
is correlated, through a canonical transformation, to another Hamiltonian system 
which lives on the standard Hilbert space of the free wave equation, the symplectic 
form not being however the usual one. 



- I/*(R 3 ) is the Hilbert space of square integrable, divergence-free, vector 
fields on R 3 . 

- M is be the projection from L§(R 3 ), the Hilbert space of square integrable 
vector fields on R 3 , onto L 2 (R 3 ). 

- (■)■)( II ' II 2 being the corresponding Hilbert norm ) denotes the scalar prod- 
ucts in L 2 (R 3 ), Lg(M 3 ), L 2 (R 3 ) and also the obvious pairing between an 
element of L|(R 3 ) and one of L 2 (R 3 ) ( the result being a vector in R 3 ). 

- Given two functions / and g in L 2 (R 3 ), f ® g is the operator in L§(R 3 ) 
defined by / ® g(A) := f(g,A). 

- H S (R 3 ), s e R, is the usual scale of Sobolev-Hilbert spaces, and iJ|(M 3 ) 
and H^(R 3 ) are defined correspondingly. 

- 9 denotes the Heaviside function. 

- I{T) denotes the compact time interval [— T, T]. 

- Lip(R 3 ;R 3 ) is the space of Lipschitz vector fields. 

- Given a measurable non negative function p, its energy E(p) is defined as 



2. Notations 
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3. The Point Limit of the Maxwell Lorentz System 

Let us consider, on the symplcctic vector space (Hl(R 3 ) x L 2 (E. 3 ) x R 6 ,f2), 

Q,((A 1 ,E 1 ,q 1 ,p 1 ), (A 2 ,E 2 ,q 2 ,p 2 )) ■= (Ai,E 2 ) - (A 2 ,E 1 ) + qi-p 2 - q 2 -pi , 

the Hamiltonian associated to the regularized Maxwell-Lorentz system in the dipolc 
approximation with an external force F = — W, i.e. 



(3.1) H r (A,E,q,p) = 2^\\E\\\ + -L ||WL||| + ^~ 



p--(p r ,A) +V{q), 



where c is the velocity of light and the square integrable density p r describes an 
extended particle with electric charge e and radius r. The corresponding Hamilton 
equations give rise to the Cauchy problem 

A r = 4irc 2 E r 



(3.2) 



E r — ~r~ AA r — o 



M (p r ,A r )p r + -f- Mp rPr 



■p r - 



1r = 75., 
p r = -VV(q) 

Mo) 

9r(0) = 90, Pr(0) 



(p r , A r ) 



A r € Hi 



E r (0)=E Q e L\ 



Po 



Let us begin with the simplest situation, i.e. V — const, and p = 0. In this case 
the fields equations decouple from the particle ones and one is led to study the 
convergence as r | of the self-adjoint operator 

Aire 2 

H r := -A H M ■ p r ® p r ■ 

m r c z 

This is a first rank perturbation of the Laplacian; its resolvent is readily calculated 
and it is given by 

(H r + Z)- 1 = (-A + Z)- 1 + r r («) _1 M • (-A + Z)- I p r ® (-A + Z)- X pr , 

where Im z ^ 0, 



T r (z) 



Aire 2 



-((-A + z)-Vr,Pr) 



and one immediately obtains that, in order to obtain a non trivial (i.e. different from 
—A) limit, the mass must be renormalized according to the classical prescription 

8ne 2 



m r := ra- 



te 2 



E{ Pr ) 



where m is the phenomenological mass. With this definition of m r , (H r + z) 1 
converges in norm to 

R m {z) := (-A + z)- 1 + TmizyHd -G Z ®G Z , 

where 



_ mc 2 V\ 
r " l(A) -~W + 6^ 



and 



1 p~Vz\x\ 

G z {x) :=-—— 
Air \x\ 



Re yfz > . 



It is not difficult to show then that R m {z) is the resolvent of a self-adjoint operator. 
In more detail one has the following result, which is no more than an adaptation 



4 



DIEGO NOJA AND ANDREA POSILICANO 



to our situation of [1, §11.1.] as regards the operator aspects and of [4] as regards 
the form ones. 

Theorem 3.1. As r J. 0, i.e as p r (x) := r p(r x), p a spherically symmetric 
probability density with bounded support, weakly converges to So, the self-adjoint 
operator H r converges in norm resolvent sense in Ll(R 3 ) to the self-adjoint operator 
H m so defined: 
1. A g D{H m ) if and only if 

3Q A eR 3 : A x :=A-^MQ A G x eH^(R 3 ), -A e p(H m ), A > , 

and the following boundary condition holds: 

1 f ( » 47re , , _ \ , mc 

hm — — / A MQ A Go ) dp r = Q A , 

rj.0 47rH J Sr \ c J e 

where S r denotes the sphere of radius r and p r is the corresponding surface measure. 
2. 

(H m + X)A := (-A + X)A\ . 

Moreover 

<Jess(H m ) = <7ac(H m ) = [0, +oo), a sc (H m ) = , 

and 

^)^|-(^) 2 }-{-Ao}, 
where — Ao has a threefold degeneration and 



X? = 2V2^n~ -M ej G Xo , 

are the corresponding normalized eigenvectors, where {ej}\ is an orthonormal basis. 
If F m is the quadratic form corresponding to H m then 

F m (A,A) + X\\A\\l = \\(-A + X)iA x \\l+(^j T m (X)\Q A \\ A>0, 

where the vector A is in the form domain D(F m ) if and only if 

3Q A eR 3 : A x :=A-—MQ A G x eHl(R 3 ). 

c 

Finally, given A G D(F m ), Q A can be explicitly computed by the formula 

3c 1 f 

Qa = limr - — - / A(x)du r . 
^ 2e r|o 47rr 2 J Sr v ; Pr 

The norm resolvent convergence of H r to H m implies, by using the explicit 
solution for abstract linear wave equations in terms of sine and cosine operator 
functions, that if 

lim \\(H r + A) 2 A r - (H m + A)5 A || 2 = 0, A > A 

rj.0 

(we used the natural distance in the energy norm between elements in the form 
domain of H r and H m ), then (see [2, corollary 2.9]) 

(3.3) lim sup \\(H r + X)i A r (t) -(H m + A)' A(i)|| 2 = 0. 

r l° \t\<T 
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Here A(t) is the solution of the Cauchy problem 

I 72 A = H m A m 



\ A{0) = A a eD{F m ), A{0) = 4nc 2 E eLt 

This gives the limit field dynamics. 

As regards the behaviour of the particle dynamics in the limit r J. 0, by [3, 
lemma 2.4], relation (3.3) implies 



lim sup 

r l° \t\<T 



This follows from (here X e Li(W) ) 



(Pr,A r (t)) - QaU) 

m r c 



. 



(p r , (H r + A) 

m r c 



-i/2 x s = e ( Pr ,(-A + \)-2X) 



2e 

37TC 



T r (x + A)- 1 ((-A + x + A)- V, X) 



((-A + A) 1 _Pr 1 Pr) dx 



m r 



(A \ ^ 1 /*°° 
—J - ^ r ro ( a; + A)- 1 (G x+A ,X) 



Vx 



and 



lim 

r|0 



lim r r (A) =T m (A) , 

r.[0 

( Pr ,(-A + X)-iX) 



rn r 



lim 

rj.0 



((-A + A)-V,/v> 



= , 

3c 2 
^87re 2 



Therefore 



lim o r 

rj.0 



lim (p r ,A r ) = Q A 

rj.0 m r c 



uniformly in time over compact intervals. 

So, to summarize, the dynamics of the system corresponding to the limit of (3.2) 
is completely specified (in the case V = const., po = 0) by solving the abstract wave 
equation c~ 2 A = —H m A and then recovering the time evolution of the particle 
position from the relation q = Qa- Since the system is linear this can be done 
explicitly (see [2, §3]). 

Let us now consider the general situation V ^ const., p ^ 0. In this case the 
fields equations and the particle ones no more decouple. However, considering the 
field equation with an assigned path 1 1— » p(t) not depending on the dynamics, one 
is led to study the Cauchy problem 

f ±A\ = -H r A r + ^f-Mp Pr 

\ A(0) = Al e #i(R 3 ), i(0)=i ei 2 (R 3 ) . 

As in the case p = one would like to show that the solution of the above system 
converges to the solution of the Cauchy problem 



$A = -Hf n A 

A(0) = A e D(F m ), 



A(o) = A ELt 
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where H m is a (not necessarily linear) operator to be identified. To this end let us 
consider the system 

f j I A r = -(H r + \)A r + ^ c Mpp r 



\ A(0) = A r e Hl(R 3 ), A(0) = A eL 

where the parameter A > A is inserted only to make H r + A, and successively 
H m + A, invertible. Its ( mild ) solution is readily found and, after an integration 
by parts, can be rewritten as 

A r (t) =cos(ct(H r + X)^)A r + sin(ct(H r + X)^)(H r + A)~^i 

Aire f* i _i 

/ coa(c(t- s)(H r + X) 2 )(H r + X) Mp(s)p r ds 

m r c J 

AlTP i 47TP l 

+ (H r + X)~ Mp(t)p r cos(ct(H r + \)2)(H r + X) Mp(0)p r . 

By [3, lemma 2.5] 

i r i 

= 



lim 

rj.0 



— (H r + Xyhdppr + -T m (X)~ 1 (H m + xyM P G x 
m r c e 



2 



and so, by thm. 3.1, if 



lim \\(H r + XyA r -(H m + XyA \\2 = , 

rj.0 

then 

lim sup \\(H r + X)^A r (t) - (H m + A)iA(t)|| 2 = , 

rj.0 | t |< T 

where 

A(t) =cos(ct(H m + Xy)A + s\n(ct(H m + Xy)(H m + A)"*i 

+ -r m (A)- 1 / cos(c(t-s)(H m + Xy)Mp(s)G x ds 
e Jo 

--T m (X)- x Mp(t)G x + ^(A)- 1 cos(ct(H m + X)i)Mp(0)G x . 
Finally, again integrating by parts, it is easily seen that 

A p (t) :=A(t) + ^r m (A)- 1 M P (t)G A 
solves the Cauchy problem 

f ± I A p = -(H m + X)A p + ±T m (X)- 1 MpG x 

\ A p (0) = A + fr„ l (A)- 1 Mp(0)G A , A p (0) =A + f r m (A)- 1 Mp(0)G A , 

and so A(t) solves the Cauchy problem 

f ^A = -(H m + X)(A+lT m (X)- 1 MpG x ) 
\ A(0) = A e D(F m ), A(Q) = A e Ll(R s ) . 

This induces us to define, on the domain 

D(H m ) := D(H m ) - -T m {\)-^MpG x , 
e 

the affine operator H m according to 

(H m + X)A := (H m + X)(a+ ^(A)" 1 M P G X ) . 



DELTA INTERACTIONS AND ELECTRODYNAMICS OF POINT PARTICLES 



7 



Alternatively (see [3, lemma 4.1]) can be defined in the following way: 
Definition 3.2. l.ig D(HP) if and only if 



A, 



A 



Aire 



MQaG\ e H 2 (R ) , 



-A e p(H m ), A > , 



and the following boundary condition holds: 



1 



lim A I 



Aire 



mc 



A MQ A G d[i r = Qa + -P ■ 



(HP l + X)A:= (—A + \)A\ . 

Note that, as it is evident, the affine operator reduces, for p = 0, to the 
linear operator H m describing the standard point interaction. 

The same considerations leading to the definition of give, coupled with a 
fixed point argument and estimates uniform in r, the following 

Theorem 3.3. (see [3, thm. 3.4]) Let V such thatVV € Lip(R 3 ;R 3 ) ; |W(x)| < 
K(l + \x\), A > A 0; and E e L^(R 3 ). Let A r e Hl(R 3 ), A € D(F m ), such that 



(3.4) 



lim \\(H r + \)'A r -(H m + \)zA \\ 2 = . 

rj.0 



Then there exists T > 0, not depending on r, such that, denoting by 

(A r , E r , q r , Pr ) e C(L(T);Hl{R 3 )) x C(I(T); L^(R 3 )) x C 2 (I{T); M 3 ) x C 1 (/(T); M 3 ) 

i/ie unique mild solution of the Cauchy problem (3.2), one has 

lim sup ||(ff r + A)'i4r(t) - (if m + A)*A(i)|| 2 = 0, 

rj.0 | t |< T 

lim sup ||£; r (t)-£;(t)|| 2 = o, 

r l° |t|<T 

lim sup \q r (t) - q(t)\ + sup \q r (t) - q(t)\ = 0, 

r l° \t\<T \t\<T 

lim sup \p r (t) -p(t)\ + sup \p r (t) -p{t)\ = , 

r l° |t|<T |t|<T 

w/iere 

(A,E,q,p) e C(I(T);D(F m )) x C(I(T);L-l(R 3 )) x C 1 (I(T);M. 3 ) x C 1 (I(T);M. 3 ) 

denotes the unique mild solution of the Cauchy problem 

A = 4nc 2 E 
E = -±W m A 

Q = Qa 
p=-VV(q) 

A(0)=A a eD{F m ), E(0) = E eL 2 ' 
<?(0) = go, p(0) =po ■ 



(3.5) 



An alternative description of the limit dynamics defined by the previous system 
is provided by the following 
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Theorem 3.4. (see [3, thm. 4.2]) Given V such that W e Lip(M 3 ;M 3 ), let 
(A,i?, g ,p)eC 1 ([0,r];D(^„))xC 1 ([0,r];L2(M 3 ))xC 2 ([0,r];R 3 )xC 1 ([0,T];M 3 ), 

A a e D(H^), E € D(F m ), be the unique strict solution of the Cauchy problem 

(3.5) . Then 

A(t)=A f (t) + ^MA s (t) , 

where Af(t) is the solution of the free wave equation with initial data Aq,Eq and 
As is the retarded potential of the source Q a 5q, i.e. 

... 1 9{ct- \x\) 
A s (t,x) = — j^j Q A {t-\x\/c) ■ 

Moreover Qa satisfies the equation 

3 C 2 3 C 3 

(3.6) Q A (t)=cVXoQ A{t) + —A f (t,0)-^p{t) . 

Let us sketch the proof of the above theorem. Given an arbitrary function Q(t), 
consider the function 

A{t, x) := A f {t, x) + ±- g(C *~ N) Q(* - \x\/c) . 

It solves the distributional equation 

1 •• 4ire 

—A = AA + MQ6 . 

c l c 

Kirchhoff formula shows that Af gives no contribution to Qa (sec [1, §3]) and so 
Q A (t) = Q(t)- Moreover, by [3, lemma 4.1] and [2, thm. 3.3] there follows 

A-7TP 

A-—-MQG X € i£(R 3 ) . 



Therefore if A(t) G D , since 



\A = AA+ —MQ5 = a(a- —MQG x ) + \—MQG x = -H^A , 
<r c \ c J c 

then the thesis will follow from unicity of the solution of (3.5). The conditions on 

Q(t) leading to A(t) e D (h^) are found as follows. By an elementary integration 



lim— ^ I (a- —MQGo 

rio Airr 2 J Sr V c 



djJL r {x) 



=A f{m + l^ Q{ t-^- Q V 



--A f (t,0)---^Q(t) 



3 c 2 rio \ r/c 
2 e 
3? 

and so A satisfies the boundary condition in Definition 3.2 if and only if Q(t) solves 
(3.6). 

Let us remark that the above theorem gives the connection with the traditional 
description: indeed the field variable satisfies a standard wave equation with a point 
source and the equation satisfied by Qa (i-e. q) is the integrated version of the ALD 
equation 

mq(t) = ^(t)-- c A f (t,0) + F(q(t)). 
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Moreover the above argument also gives the converse statement: 
the solution of the distributional Cauchy problem 

±A = AA + ^Mq6 
q(t) = c^r Q q(t) + ^ A f {t, 0) - ffj p{t) 
p = -VV(q) 

A(0) = A a e D(H%), A(0) = A7rc 2 E e D(F m ), 
q(0) = q , q(0) = Q Ao , p(0) = p 

solves (3.5). Also note that the equivalence between the two descriptions holds true 
if and only if the initial data for the field are chosen coherently with the particle's 
ones (i.e. q a = Qa and % = Qa„)- 

Let us now come to the Hamiltonian character of system (3.5). By definition 
3.2, and by D(-ffPJ C D(F rn ), one can check that 

(HP l A 1 ,A 2 )=F m (A 1 ,A 2 )+47rp-Q A2 . 

Therefore equations (3.5) are nothing but the Hamilton equations corresponding to 
the (degenerate) Hamiltonian 



H m (A, E, q,p):= 2irc 2 \\E\\ 2 + i- F m (A, A) + p ■ Q 

this is defined on the symplectic vector space (D(F m ) x L 2 
one has the following convergence result: 



V{q) ; 

R 6 ,fi). Moreover 



p--( Pr ,A) +V(q) 



(q,p) e M. 6 , and let 



Theorem 3.5. (see [3, thm.4.5]) Let 

H r {A, E, q,p) = 2kc 2 \\E\\ 2 + i_ ||VA||* + -L 

be the Hamiltonian giving the equations (3.2), let E £ L 2 
A r G Hl(R 3 ), A e D(F m ) satisfy the condition (34). Then 

lim H r (A r ,E,q,p) =H m (A,E,q,p) . 

rj.0 

4. The Non Runaway Dynamics 

The negative eigenvalue in the spectrum of the operator H m gives rise to un- 
stable behaviour which corresponds, in classical electron theory, to the presence of 
the so called "runaway solutions". In this section we briefly describe, in the free 
case and, for simplicity of presentation, vanishing particle momentum, the reduced 
dynamics on the stable manifold. 

Given any vector subspace V C L 2 (M. 3 ), we define the corresponding "non 
runaway" subspace [V] nr by 

[V] nr :={AeV : (A,X$)=0, j = 1,2,3}, 

X° being the eigenvectors corresponding to — A (see thm. 3.1). Observe that, if 
A = A x + ^ MQ A Gx e D(F m ), then one has 

Aire 
c 

-c 



(A, X?) ={Ax,X$) + — (MQ A Gx,X°) 



=2V2 



irm- 



(AlG Xo ) + ^Q 3 A (Gx,Gx ) 
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Therefore A e [D(F m )] nT if and only if 

3c (G Xo ,A x ) 
^ A 8ne (G\, G\ ) ' 
so that any A e [D(F m )] nr is uni vocally determined by its regular part A\. This 
implies, since (Ga ,Ga ) _1 = 8ny^Xo, that the map 

$ : Hl(R 3 ) - [£>(F m )] nr , $A := A + MP A G X(] , P A := — 127t\/Ao (G Ao , A) , 
is bijective, with inverse given by 

fc" 1 : [D(F m )] nr - Hl(M. 3 ), Q-'A := A MQ A G Xo = A\ . 

c 



Lemma 4.1. $ is a continuous bijection between H I (R 3 ) and [£>(F m )] nr . More- 



over 



(4.1) $- 1 ([^(i/ m )]„r) = {4 € ff?(R 3 ) : A(0) = 0} . 

Proof. If A > A we can use \\A\\\ := F^(A,A) as a norm on [D(F m )] ni . Then 



= \\A + MP A (G Xo - G\) + MP A G X \\ 

2 



A 



=||(-A + A)* (A + MP A (G Xa - G x ))\\l + - T m (X)\P A \ 2 

<2(||(-A + X)iA\\l + ||(-A + A)^MF A (G, - G A )|| 2 ) + | T m (X)\P A \ 2 

<2||(-A + A)^||2 + C1 |P A | 2 < 2\\(-A + X)iA\\ 2 2 + c 2 \\A\\ 2 2 

<c 3 ||(-A + A)^||i. 

Since, for any A ^ A , A = A x + T m (X)- 1 MA x (0)G x G [£>(ff TO )] nr if and only if 
-4a (0) = T m (X)P Ax , one has 

($" 1 A)(0) =A A (0) + (M P A (G x - G Ao ))(0) 

=A A (0)-r m (A)P j4A = . 
Then is continuous by bijectivity and by the open mapping theorem. □ 

Obviously the map $ can be extended to the whole £^(R 3 ), giving rise the orthog- 
onal projection onto [L^(IR 3 )] nr . 

Let us now consider the symplcctic space ([D(F m )] nr x [Z^(R 3 )] nr , ri ) , where £1 
denotes the canonical symplectic form induced by (•, •), i.e. 

n ((A 1 ,E 1 ),(A 2 ,E 2 )) = (A U E 2 ) - (A 2 ,E 1 ) . 

On ([D(F m )] nl x [i^(R 3 )] nr , fio) we have the non negative Hamiltonian 

H+(A,E) = 2nc 2 \\E\\ 2 + ±-F m (A,A) 

2 \ 



^\\E\\l + ± HVA Ao || 2 + ^(^) Q,\ 2 



with corresponding Hamiltonian vector field 



X n +(A,E) = (^nc 2 E,-^-H m A^j , 
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defined on the domain [D(H m )] ni x [D(F m )] nr . If we pull-back Slo and to 
Hl(WL 3 ) x Ll(R 3 ) by using the map 

* := $ x $ : Hl(R 3 ) x L?(R 3 ) - [D(P ro )] nr x [L2(R 3 )] nr , 

we obtain the following 

Theorem 4.2. 1. (i^(R 3 ) x Ll(M. 3 ),Sl ni ) is a symplectic space, where the 
(weakly) nondegenerate symplectic form £l nr is given by 

n nr ((Ai, Ex), (A 2 ,E 2 )) := **n ((Ai, Ex), (A 2 ,E 2 )) 
=n ((Ax,Ex),(A 2 ,E 2 ))- 1 (P Al -P B2 -P A2 -PeJ . 



Defining 

X Hni : fl?(R 3 ) x ^(M 3 ) -> Ht(R 3 ) x 



where 



1h„ (A, £) := ^WS, — AA + - ^X MA(0)G Xo 
one has, for any (A, E) G H?(R 3 ) x iT^R 3 ), 

1 o nr (x„ nr (A e), (A, £)) = n ni (A, E) , 

H m :Hl(R 3 ) x Ll(R 3 ) — > R 
W nr (A,P) :=**W+(A,.E) 

1. e. X-H nr is ifce (unique) Hamiltonian vector field corresponding to H. nl . 
3. The Hamiltonian vector fields X n + and X-^ nr are VP -correlated, i.e. 

(4.2) X w + o * = * o X Hnr , 
and 

(4.3) E/ m (t)otf = tto£7 nr (t) , 

w/iere U m (t) and U nl (t), t £ R, denote the one parameter groups of canonical 
transformation of ([D(F m )] m x [L^(R 3 )] nr , Oo) a™d (#*(R 3 ) x £*(R 3 ), O nr ) given 
6y i/ie flows of X H + and Xfi ni respectively. Moreover, ifUf(t), denotes the 

flow, on Hl(R 3 ) x L^(R 3 ), given by solving the free wave equation, then 

( 4 - 4 ) U m (t)\ [r>(i? m )] nr x[r>(F m )] nt = * ° E/>(t) ° • 

Proof. 1. A simple calculation shows that 

n ({$Ax,$Ex),(<S>A2,$E 2 )) = ($Ax,$E 2 ) - ($A 2 ,$Ex) 

=n Q ((Ax,Ex),(A 2 ,E 2 ))- 1 (P^-P^-P^-PeJ , 

1Z7TV Ao 

2. Since 

= ||(-A + Ao)'A||l - \ \\$A\\l 
=||(-A + Ao)U||2-A (A $A) 

HIv^ + ^IPaI 2 , 
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\\*E\\l = (E,*E) = \\E\\l-—± 7 =\P E \* , 

1Z7TV Ao 

posing 4irH nT A = —AA — &n\/\QMA(Qi)G\ , we need to verify the relation 
-(AA, A) + ^ \P A \ 2 = MH ni A, A) - ^-L= P 47lHnrA ■ P A . 

Since 



and 



one has 



P AA = - 127rv/Ao(-(G Ao , (-A + X )A) + A (G Ao , A)) 
=12ny/j^A(0) + X Q P A , 

Pma { o)g Xo = -12ttv^| —L^= A(0) = -A(0) , 



47r(ff nr A, A) - 1 P 4 nH niA ■ P A 
= -{AA,A)- l -A{Q)-P A 



+^-L= (l2^v^A(0) + \ P A ) -Pa- I A(0) ■ P A 



-^0 m ,2 



= - { AA,A) + ^\P A 

3. Xft nr generates a group of continuous (w.r.t. the Hilbert norm on Hi 
L^(R 3 )) linear transformations since we can write 



47rH nr + A = (1 - 67t\/A Ga ® G Ao ) o (-A + A ) . 

Moreover such transformations are symplectic, since X-h^ is Hamiltonian. Formula 
(4.2) follows from the definitions of fi nr and Ti nri and (4.3) follows from (4.2) by the 
unicity of generators. Finally, (4.4) follows from (4.1) and the definition of X~n . 

□ 
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